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Non-Perturbative Renormalization of Lattice QCD. 

G.C. Rossi a * 

a Dip. di Fisica, Universita di Roma "Tor Vergata" and I.N.F.N., Sezione di Roma II, 
Via della Ricerca Scientifica 1, 1-00133 Roma - Italy 

In this talk I will discuss a number of approaches designed to deal with the problem of setting up a fully 
non-perturbative renormalization procedure in lattice QCD. Methods based on Ward-Takahashi identities on 
hadronic states, on imposing chiral selection rules on amplitudes with external quark/gluon legs, on the use 
of the Schrodinger functional and on "heating and cooling" Monte Carlo steps are reviewed. I conclude with 
some remarks on the possibility of defining next order terms (higher twists, "condensates", ...) in short distance 
expansions. 



1. Introduction 

Renormalization is a necessary step to bring 
numbers extracted from Monte Carlo simulations 
in contact with actual physical data. Its role is 
threefold: 

• to allow the construction of finite operators; 

• to recover (modulo possible anomalies) chiral 
symmetry, broken by lattice regularization; 

• to connect the high momentum perturbative 
regime of QCD with the low momentum non- 
perturbative region of the theory through the run- 
ning of renormalized quantities. 

Perturbative calculations of Renormalization 
Constants and Mixing Coefficients (RC&MC's) 
can only be of limited value in numerical simula- 
tions, because 

• uncontrollable Non-Perturbative (NP) contri- 
butions may affect dimensional MC's, 

• only the lowest order terms of perturbative 
expansions can be calculated, 

• field theoretical perturbative series are (most 
probably) only asymptotic. 

Methods to compute RC&MC's beyond Per- 
turbation Theory (PT) are thus highly desirable, 
if not for the fact that, in a strict sense, in a 
fully NP approach of QCD use of PT should 
nowhere be made. In the recent past, and es- 
pecially in the last year, a number of significant 
progresses have been made in the direction of de- 
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veloping efficient strategies to compute RC&MC's 
in a NP way, i.e. directly from appropriately de- 
signed Monte Carlo simulations. The plan of the 
talk is the following. In Sect. 2 I start by recall- 
ing how Ward-Takahashi identities (WTI's) on 
hadronic states can be used i) to define partially 
conserved SU(Nf)L <8> SU(Nf)x currents, obey- 
ing Current Algebra (CA) and ii) to construct 
finite composite operators with well defined chi- 
ral transformation properties. In Sect. 3 I dis- 
cuss the NP determination of RC&MC's, based 
on imposing chiral selection rules and renormal- 
ization conditions on Green functions with (am- 
putated) quark/gluon external legs. This ap- 
proach amounts to require the validity of WTI's 
on quark/gluon states. In Sect. 4 the use of the 
Schrodinger Functional (SF) for the determina- 
tion of the 0(a) improved form of the fermionic 
action and of quark bilinear RC&MC's is illus- 
trated. For completeness I will briefly recall in 
Sect. 5 the "heating-cooling" method designed to 
extract from pure gauge Monte Carlo simulations 
the RC of the topological charge density and the 
divergent subtraction needed to arrive at a finite 
definition of the topological susceptibility. Fi- 
nally I conclude in Sect. 6 with a few observations 
on the problem of giving a rigorous NP definition 
of large scale effective expansions beyond leading 
order and on the related question of the possibil- 
ity of evaluating the ensuing higher dimensional 
corrections. 
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2. WTI's on hadronic states 

For the purpose of this talk lattice QCD should 
be considered as a regularized version of the the- 
ory described by the (euclidean) continuum La- 
grangian 



C-QCD = £ym + C-f — ^^{F^vF^v)- 



(1) 



As is well known, to avoid fermion species dou- 
bling a dimension 5 (irrelevant) operator, the so- 
called Wilson term [y , must be added to the naive 
lattice discretization of ([[]) , leading to the follow- 
ing expression for the fermionic part of the action 



^E^(zW*)(r-7M)V>(z + |u) 



The Wilson term (r J cfixipD 2 ^ in continuum 
notations) explicitly breaks global SU(Nf)L <g> 
SU(Nf)n chiral symmetry. Chiral symmetry is 
an exact property of the continuum Lagrangian, 
which is only softly broken by quark mass terms. 
Actually the Lagrangian ([[]) possesses two further 
U(l) symmetries: an exact U(l)v vector sym- 
metry associated to the conservation of fermion 
number, and an axial U{1)a symmetry, which is 
explicitly broken at quantum level by the color 
anomaly and is responsible for the singlet/non- 
singlet pseudoscalar mass splitting. The same 
symmetry pattern is seen to emerge from the 
a — y limit of the regularized theory. 

2.1. Vector and axial currents 

It has been shown, in fact, in refs. || || that in 
the continuum limit chiral symmetry can be fully 
recovered. Partially conserved vector and axial 
currents, V* and A^, / = 1, —,Nj — 1, satisfying 
CA can be defined. In lattice QCD finite RC's 
relate, in the limit a — > 0, operators obeying CA 
to their bare expressions 



Vj = ZyVf , 



4 = Z A A{ 



(2) 



The RC's Zy and Za can be determined non- 
perturbatively precisely by requiring the valid- 
ity of the non-linear constraints coming from 



CA g| (J] [|. Once Z A has been determined, im- 
posing PCAC in the (naive) continuum form fixes 
the linearly divergent mass subtraction term, M, 
which allows to make the pseudoscalar operator, 
^ 5 {X f , M - M}ip, finite. 

2.2. Other composite operators 

The explicit breaking of chirality, induced by 
the Wilson term, also spoils the "nominal" chi- 
ral properties of composite lattice operators. Let 
0[ n ] = {0[ n ]} De a basis of operators which 
naively (i.e. at the tree level) transform accord- 
ing to the irreducible representation, [n], of the 
chiral group 



i Sat - [r M> 



(3) 



where I have denoted by SO 1 ^ /Sat the variation 

of O) , under an infinitesimal global chiral rota- 
[n] o 

tion and by r<, the f th generator of the chiral 
group in the representation [to]. 

Explicit lattice perturbative calculations § 
show that, as expected, radiative corrections 
induce mixing among operators with different 
(nominal) chiral properties. The existence of par- 
tially conserved lattice vector and axial currents 
guarantees, however, that, given the set of bare 
operators O^j , it will be possible to determine 
the coefficients, cV /n , of the linear combination 



01 



01 



E 



c n° 



(4) 



in such a way that, up to terms of 0(a), 0| n j 
will obey the set of WTI's appropriate for an op- 
erator belonging to the representation [to]. The 
resulting operator will automatically be multi- 
plicatively renormalizable. In the indices [to'] 
and j run over all operators with dimensions 
equal or smaller than and belonging to all 
possible representations of the chiral group with 



2 Note that the expression of the bare currents is not 
uniquely determined, but can be modified by the addi- 
tion of higher dimensional operators, vanishing in the limit 
a —y 0, with the same conserved quantum numbers as the 
currents. In doing so the RC's Zy and Za will have to be 
accordingly modified. A similar observation holds for any 
lattice operator. 



3 



in the limit 



the only constraint of having the same conserved 
quantum numbers as . The mixing coeffi- 
cients of operators with the same dimension as 
0[n] are dimensionless and finite, while the coef- 
ficients of the lower dimensional ones will diverge 

ij -(dimof„]-dimof 
as c, ~ a i' 

[71,71 \ 

In this talk for lack of space I will limit my 
considerations to the massless case. Since flavor 
vector symmetry is unbroken by the lattice reg- 
ularization, to determine the c's it will suffice to 
look at the axial WTI's. One way to proceed is to 
impose that the renormalized, integrated, lattice 
axial WTI 



x{(x))Ol n] (0))\h 2 ) 



-i(hi 



l *Q[n](°) 
I 8^7 



h 2 ) 



(•5) 



taken between on-shell mesonic states, \h\) and 
\h%) [], has the form expected for an operator be- 
longing to the representation [n], i.e. that 

EMTixfapUO)) \h 2 ) \h 2 ) 



= (rf n] ) tf (^|Of w] (Q)l^) 

By varying the external states, one can write a 
sufficiently large number of (linear) equations and 
fix the c's uniquely. will be finally made 

finite multiplying it by an overall RC defined, for 
instance, by the condition 

(fti|0[„](/i) |/i 2 >| M 2 = {hi\Z(aii)d [n] (a)\h 2 )\^ = 
= continuum normalization 

In practice this procedure can only be employed 
for bilinear quark operators (such as currents, 
scalar and pseudoscalar densities,...). For more 
complicated operators, like the four-quark opera- 
tors describing the Effective Non-Leptonic Weak 
Hamiltonian (ENLWH), this approach would re- 
quire measuring with practically unattainable 
precision a much too large number of matrix ele- 
ments In the following sections I will describe 

3 As chiral symmetry is spontaneously broken, the axial 
rotations of the Of ll and Oh 2 operators, which create the 
mesons hi and /12 from the vacuum, do not contribute to 
on-shell matrix elements. 

4 See, however, the method proposed in g|, based on the 
idea of fixing the power divergent MC's from the relations 



alternative strategies, recently proposed to over- 
come this kind of difficulties. 

3. WTI's on quark and gluon states 

The simple idea of is to compute RC&MC's 
mimicking in Monte Carlo simulations the 
straightforward procedure employed in contin- 
uum PT. To explain the method let me consider 
the physically interesting case of the AS = 2 four- 
quark operator O as=2 = {sj^d)(sj^d), whose 
matrix elements control the K —> 7r7r ^7=3/2 de- 
cay amplitudes 0] ||]. In lattice QCD this oper- 
ator mixes with 4 other operators of dimension 6 
with dimensionless coefficients. The finite opera- 
tor, O as=2 , has the general expression 



qAS=2 _ ^AS=2^qAS=2 



AS=2qAS=2~ 



(6) 



The proposal of is to extract from Monte Carlo 
data the (sd\Of° \ds), i = 0,1,. ..,4 matrix el- 
ements and to fix the coefficients Z AS=2 , by re- 
quiring that, at large p\ = ^ 2 (the pkS are the 
momenta of the external legs), all chirality vi- 
olating form factors of the renormalized opera- 
tor, O as=2 , should be identically zero. In prac- 
tice this means that in the interacting theory 
qAS=2 mus t precisely match, at large /j, 2 , the fla- 
vor, color, spin,... structure of the bare operator, 
O as=2 . The overall RC can be successively fixed 
by using, for instance, the same renormalization 
condition employed in the continuum. Techni- 
cally the whole procedure is carried over by first 
defining the matrix (i, j — 0, 1, ...,4) 



D ii ( M )=Tr(P i Af s=2 (M9;M)) 



(7) 



where the amplitude A AS ~ 2 (qqqq; fj.) is the four- 
legs amputated matrix element of the operator 

O: 



AS=2 



A 



AS=2 



with each leg taken at momentum p = fx 
(MM) = {q{p)q{p)0 AS = 2 q(p)q{p))\^ 



The P's are orthogonal projectors, Tr(PiPj) = 
Sij, satisfying Tr(PiA AS=2 (qqqq)\ tree ) = 

coming from the low energy theorems of the chiral sym- 
metry (Soft Pions Theorems - SPT's), while computing in 
PT the overall RC and all the other dimensionless MC's. 
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MC's are determined from the constraints [] 
Tr(PiA AS = 2 ) =0,i= 1,...,4 

A A5=2 = A AS=2 + £4^ Z AS=2 A AS=2 



IS 



(8) 



by solving the non-homogeneous set of linear 
equations 

4 

Y / Z^ 2 D ji = -D 0i ,i=l,...,4 

3 = 1 

The overall RC is evaluated from the condition 



zAS=2 



VTr(P A 



AS=2\ 



(9) 



where Z g is the quark wave function RC. 



In 



four Z q 1 ' 2 factors appear, because the A's are 
four-leg amputated amplitudes. In Fig. 2 of Q 
the results of a number of Monte Carlo measure- 
ments of the (K \O AS=2 \K) amplitude are shown 
as a function of m 2 K @ . It is clearly seen that the 
long sought chiral behaviour of (Ko\O as=2 \K) is 
correctly reproduced when the NP renormaliza- 
tion procedure described in this section is em- 
ployed, in conjunction with the SW-Clover-Leaf 
improved fermionic action || . Similarly good 
results have also been reported by the JLQCD 
collaboration in the case of the standard Wilson 
action fll|| . 

A possible difficulty with the idea of extract- 
ing the NP values of RC&MC's from quark/gluon 
matrix elements is that a lattice gauge fixing is 
required in the simulations, leaving behind the 
problem of Gribov ambiguities. 

The next obvious step within this approach is 
to go over to the much more complicated and in- 
teresting case of the AI = 1/2 ENLWH, where 
the problem of the subtraction of lower dimen- 
sional operators with power divergent MC's has 
up to now forbidden a reliable computation of the 
K — * 7T7r| A/=1 / 2 amplitude |l2| |13|. The two op- 
erators relevant in this case belong to the [8, 1] 
representation of the chiral group. They are usu- 
ally indicated by and their bare expression 



5 In the equations below for notational simplicity we drop 
the arguments of the amplitudes A whenever unnecessary. 



o, 



(-) 



+2{s 1 L d){d 1 ^d) + 2(srfd)(srfs)}- 

As we repeatedly said, the explicit breaking of 
chiral symmetry due to the presence of the Wil- 
son term in the action, induces the mixing of Oq 
with operators belonging to chiral representations 
other than the [8,1]. Taking into account the 
symmetries left unbroken by the lattice regular- 
ization, it is easily seen that the renormalized, 
finite lattice operators, which in the fully inter- 
acting theory will transform as an [8, 1] repre- 
sentation, must have the expression 



0(±) =Z (±) [0 (±) +E 4 =iZ (±)^(±)_ 

+Z^ ± 'sa /J . I/ F^ v d + z\ sa^F^d-A 



(10) 



sd - 



tP{±: 



sj 5 d] 



Let us in turn examine the various terms in (|l0|). 

• Dimension 6 operators 

The spin and color structure of the operators 
of dimension 6 contributing here is obviously the 
same as the one of the AS — 2 case discussed 
before. Only the flavor structure is different. 

• Dimension 5 operators 

If the GIM mechanism is operative, the coeffi- 
cients Z§ are actually finite, because the po- 
tential 1/a divergence is replaced by a m c — m u 
factor. As for Z 5 P(±) , GIM and CPS symme- 
try |Q (CPS = CP x symmetry under s — » d ex- 
change) make it finite and vanishing in the limit of 
exact vector flavor symmetry (exact SU(Nf)v)- 

• Dimension 3 operators 

The GIM mechanism softens the divergence of 
.Z3 , reducing it from 1/a 3 to (m c — m u )/a 2 . As 
before, CPS symmetry makes Z 3 vanishing in 
the limit of exact SU(Nf)v- 

An important observation at this point is 
that the Maiani- Testa no-go theorem jl4| (which 
states that in euclidean region essentially only 
matrix elements of operators between one-particle 
states can be extracted from Monte Carlo data) 
forces us to limit our considerations to the 
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{K\b^ \ir) matrix elements, leaving the recon- 
struction of the K — ► 7T7T amplitude to the use of 
the SPT's 



10; 



8,11 



i2a 2 



(n + (p)\0^ 1 \K+(q)) = -2a 2 



/ + _i / -,AS=l| 



'[8,1] l A / — 44 Ql — F| — 

valid for an [8, 1] AS = 1 chiral operator. From 
the above equations it is in fact immediately seen 
that the K — > nn amplitude can be obtained from 
the slope of (tt+^IO^ 1 as a function 

of (p ■ q). Notice that, proceeding in this way, 
only the MC's of the parity-conserving operators 
contributing to (10) will be necessary. 

Let me now briefly illustrate the strategy for 
the construction of the AI = 1/2 ENLWH (see 
also Q). The method we propose to arrive at a 
NP evaluation of the finite and power-divergent 
MC's in jl0| ) is a generalization of the approach 
described before to deal with the case of O as=2 . 

To avoid uncontrollable numerical instabilities 
when the simultaneous computation of finite and 
infinite (in the limit a — > 0) MC's from a sin- 
gle (large) set of linear equations is attempted, 
it is more convenient to separate the subtraction 
of power divergences from the subsequent finite 
mixing. We thus subdivide the whole procedure 
in two steps. 

• One first defines the "intermediate" (power 
divergent) mixing constants, C§ 
by requiring 



^ sd) and Cf d)(±) 



Tr(P s - d Ag ) feM)) = 0, i = 0,l, 
Tr(P s - d A 5 (g 9 ;/i))=0 



(11) 



where P S d is the projector over the color and spin 
structure of the operator sd and the amplitudes 
Ag^(gg; /i) and h.§(qq\ /x) are respectively the am- 
putated q, q matrix elements of the operators 



O 

o 5 



(±) 

6,i 



o 



(±) 



a 



sa 



(sd)(±) 

(M) - 



sd, i = 0, 1, ...,4 



^F^d+cr'sd 



• The second step consists in determining the 
finite MC's, and by imposing the (non- 
homogeneous) set of linear conditions 

T*(P 6 A<g(g ff g;/x)) 

r(±) 



o 



Tr(P i A^(qqqq ]fJ ,))=0, i = 1, 



where the amplitudes A^ are the appropriately 
amputated matrix elements of the (by now only 
logarithmically divergent) parity-conserving part 
of the operator dlQ), i.e. of the operator 



U PC — 



PC 



A feasibility study of the procedure outlined 
above is under way EEI . 

4. Schrodinger Functional 

The functional integral with Dirichlet bound- 
ary conditions along the time direction 



/C($ 2 ,*i;T) = 



$ 1 (x)=$(a,0) 



(12) 



yields the amplitude to find the field configura- 
tion $2 = $2(%) at time t = T, if the field con- 
figuration was $i = &i(x) at time t = 0. It thus 
represents the matrix element of the (Euclidean) 
transfer matrix, e~ HT , between the Schrodinger 
states |$ 2 ) and |$i): 



/C($ 2 ,$i;T) = ($ 2 |e- HT |$i) 
More generally, given the functionals 
Qfe(xfe) = Qk[®{2k,tk),8/6®(x k ,t k )], k = l, 
one can compute the expectation values 

(Ql{xi)Q 2 (x2)-.Qn(x n )) = 

= f* 2 V^]e- s ^Q 1 (x 1 )Q 2 (x 2 )...Q n (x n ) 
= (®2\e- HT nLi e Ht *Q k (x k )e- Ht * \^) 



(13) 



In QCD one must take $ 2 = {A 2 , ip^ ' i '} and 



$i = iVi }) where the superscripts 

(±) appended to the "tp,ip fields are there to re- 
member us that only half of the fermionic compo- 
nents can be assigned at each time boundary []l6f . 

The YM and the QCD SF's have been formally 
constructed in the continuum, in the temporal 
gauge (Aq — 0), in [jl7| and in 18 respectively. 
In Q the superscripts (±) were taken with ref- 
erence to the positive and negative frequency de- 
composition of the fermionic fields, as this is the 
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most suitable choice when working in the tem- 
poral gauge. The corresponding constructions on 
the lattice, where no gauge fixing is necessary, 
were carried out in (l9) and in pOfl . In the QCD 
case Dirichlet boundary conditions on the spin 
projections 



(±) 



1±7Q 
2 



1p, = 1p 



l±7o 
2 



were imposed pOfl . Studies of the renormalization 
properties of the lattice QCD SF can be found 
in |lj and in ||. 

In absence of fermions, the SF formalism was 
implemented in lattice simulations to study the 
running of the gauge coupling. A renormalized 
coupling, a s , as a function of q = 1/L (L d = lat- 
tice volume) is defined by looking at the response 
of K. under small changes of an externally applied 
color-magnetic field |23| ] . The running of a s with 
1 / L is extracted in a recursive way, in a two step 
procedure 0]. 

• Starting with a given (small) value of g$, the 
variation of a s , Actf = a s (L) — a s {L\ — sL) 1 
ensuing from an increase of the number of lattice 
points by a factor s d (L = Na — ► L% = sNa), is 
computed and its continuum limit is extrapolated 
from data taken at decreasing values of a/L with 
fixed a s (L). 

• Then the physical value of the lattice spacing 
is increased, at fixed renormalized coupling (by 
adjusting the bare coupling g^), in order to have, 
with a number of points equal to the initial one, a 
renormalized coupling precisely equal to a s (Li), 
to be in position to start a second iteration. 

At this point the whole procedure is repeated, 

(2) 

by measuring the variation Acts = a s (Li) — 
cts{L2 = sLi) for an increase of the number of 
lattice points again by a factor s d . At each step 
the limit a/L — » of the variations Acvg is ap- 
propriately taken. 

The results obtained with this method are 
very accurate (few % errors) and compare very 
nicely |^5| with the results of 26 1, derived em- 
ploying a similar iterative procedure, but using 
for the renormalized gauge coupling a definition 
given in terms of ratios of twisted Polyakov loops. 
For a measure of the running of a s defined di- 
rectly from the three-gluon vertex, see |27) . 



The introduction of fermions in the SF opens 
the way to a wealth of very accurate NP mea- 
surements of RC&MC's which ultimately will al- 
low the (on-shell) construction of the fully 0(a) 
improved QCD lattice theory. The use of the SF 
approach has the further noticeable advantages 
of providing a formalism which allows to perform 
Monte Carlo simulations 

• directly at the chiral point, Mo = M cr , 

• with no need of fixing the gauge. 

Up to now the method has been used to con- 
struct (in the quenched approximation) the 0(a) 
improved QCD action and the correspondingly 
improved bilinear quark operators. However, the 
approach seems to be sufficiently general to be ca- 
pable to encompass the much more complicated 
and interesting case of the four-quark operators. 

The general idea is to start with the SW 
nearest-neighbor improved fermion action 



'-c V 



where S w is the standard Wilson action (from 
now on r = 1) and is any lattice discretiza- 
tion of the gauge field strength |2gj] , and to de- 
termine the value of c sw and of all the neces- 
sary RC&MC's, by requiring that on-shell chiral 
WTFs involving quark bilinear do not have 0(a) 
corrections. The "tree level" improvement of the 
Wilson action (eq. ( |l4| ) with c sw — 1) kills in on- 
shell Green functions all terms that in PT are of 
0(a (<7q log a) n ) (i.e. the terms that in the contin- 
uum limit, <7g ~ 1/loga, are effectively of order 
a), leaving uncancelled 0(a) subleading logarith- 
mic corrections JlC) ] . 

To illustrate the SF method let me describe the 
strategy for the determination of quantities like 
c sw , Za, Zy or the critical mass, M cr , which are 
only functions of g$ . More subtle is the problem 
of computing the renormalized coupling constant, 
g\, or the renormalized quark mass, mjt, or cer- 
tain MC's, because it turns out that the general 
mass-independent renormalization scheme, con- 
sistent with 0(a) improvement, requires a rescal- 
ing of the bare parameters by mass dependent 
factors [§§. 

The crucial observation which makes the whole 
program actually feasible is that, for the sake of 
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computing the functional dependence of c sw , Z A , 
Zy, or M cri on g%, for small g$, all the necessary 
numerical simulations can be quite happily per- 
formed deep in the perturbative region. This can 
be realized by working in a not too large volume 
in order to keep lattice momenta, p n = tt-ti, much 
larger than the typical mass scale of the theory. 
Indicatively, if N is the number of points per lat- 
tice side, one should have 

2tt/L > Aqcd i-e. N < exp(l/2& So) 

At the same time to keep under control discretiza- 
tion effects, which scale like powers of a/L (times 
possible logarithmic factors), one must also re- 
quire 

L/a = N ~ large 

With the choices g% < 1 and N > 10 both condi- 
tions are well satisfied and the physical volume of 
the box (determined, for instance, as explained 
in (23)) turns out to be rather small (L ~ 0.5 
fm). An important consequence of these choices 
is that simulations can be performed directly at 
the critical mass, Mq = M cr , as for small <?q the 
lowest eigenvalue of the Dirac operator is 0(1/L) 
and not O(M - M cr ). 

To illustrate the method in a concrete way, 
let me discuss the determination of c sw {g^) and 
M cr .(<?g). One starts from the lattice PCAC equa- 
tion 

(dpAfaxoWf) = 2m(P f (x,x )O f )+O(a)(15) 
where {x ^ Support of O*) 

, - T* , - T* 

is an operator which creates a pseudoscalar 
state at t = 0. Explicitly one can take 



of 



d 6 y / d 6 z — — 75- — ^— — 



In (|lq ) it is understood that the fermionic bound- 
ary values are set to zero after the action of the 
grassmannian functional derivatives. A plot of 

R A = (a„A^S,x )Of)/(pf(S,x )Of) 




Figure 1. c sw for large and small g\. Compari- 
son with 1-loop TIPT (crosses) and PT (dots) is 
shown. The solid line is a numerical interpolation 
of the data. 



as a function of the current insertion time, Xo/a, 
shows strong violations of chirality 29 1 , when the 
standard Wilson action is employed. These ef- 
fects are to be attributed to the (large) 0(a) cor- 
rections affecting eq. (|l5| ) . They can be eliminted 
from all Green functions at unequal points, fol- 
lowing the steps described below 

• add the SW-term, \c sw a 5 ^ a ^^P^, to 
the standard Wilson action, 

• construct renormalized improved quark oper- 
ators, O q , by adding to their bare expression, O q , 

(i) 

a suitable linear combination of operators, Oq , 

(i) 

with dimensions up to dimOq = dimOq + 1 , 

• fix, at any given value of <?q, c sw and the 
corresponding MC's by requiring 0(a) corrections 
to be absent from WTI's. 

The renormalized, 0(a)-improved expressions 
of the axial current and of the pseudoscalar den- 



sity are {m q = M — M cr ) 

Af* = Z A [(1 + am q b A )Al + ac A d^Pf] 
P fI = Zp(l + am q b P )pf 



(16) 



where at tree level, beside c 

b tree = tfree = ^ Jree = g |j 

<7o the improved ratio 



= 1, one has 
At any given 



Ra = {d»Af (x,x )Ot)/{P aJ (x,x )Of) 

is measured for different values of xq/o. and for 
two different (zero and non-zero) boundary gauge 
field configurations. Actually to the order one is 
working (terms 0(a 2 ) are neglected), the behav- 
ior of R A as a function of xq /a is only sensitive 
to c sw and ca- In fact one can write 



Z A l+am q b A ((c%Al+ac A d 2 pf)Qt) 



R 1 

A Zp l+am q bp 



0(a 2 ) 



The nice results obtained for c sw {g Q l ) and ca(<7o) 
by imposing that R A is independent of i) xo/a, 
ii) the particular form of the operator iii) the 
boundary values of the gauge fields, are shown in 
Figs. 1 and 2. While ca is negligible to all practi- 
cal extents, c sw appears to be a rapidly increas- 
ing function of g 2 . Data for c sw are larger than 
expected from simple "tadpole improvement" ar- 
guments [M. Once c sw and ca have been fixed, 
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Figure 2. ca as a function of g$. The dotted line 
is 1-loop PT. 



by varying Mo, one can determine M cr as the 
value of M for which Ra (and of course also R A ) 
vanishes. Data for n cr , measured at large and 
small (?oj are shown in Fig. ||, together with the 
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Figure 3. n cr for large and small g^. Compari- 
son with 1-loop TIPT (crosses) and PT (dots) is 
shown. 



results from 1-loop Tadpole Improved PT (TIPT) 
and from 1-loop straight PT, respectively. It is 
clearly seen that the non-monotonic behavior of 
K cr is in striking contrast with tadpole improve- 
ment expectations. This is not surprising in view 
of the fact that a proper definition of M cr re- 
quires a linearly divergent subtraction. Absence 
of O(a) corrections in WTI's is a strong require- 
ment which also allows an accurate evaluation of 
the current RC's. In Figs. [|and || we show the re- 
sults for Za and Zy, as functions of g^, together 
with the predictions of 1-loop TIPT and PT. In 
both cases numbers from TIPT interpolate quite 
nicely Monte Carlo results. Perhaps the lesson 
we can draw from the different ability of TIPT in 
reproducing data directly extracted from simula- 
tions is that TIPT can be really useful only for 
estimating finite RC&MC's. 

We would like to conclude this section by stress- 
ing again that the use of the SF formalism in 
Monte Carlo simulations seems to be extremely 
promising and has already provided us with a 
lot of very accurate results. It would really be 
a great achievement if one could extend this ap- 
proach to the construction of finite, renormalized 
four-quark operators. 
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5. NP renormalization by "heating - cool- 
ing" steps 

For completeness, in this section I wish to dis- 
cuss a completely different NP renormalization 
method based on "heating-cooling" Monte Carlo 
steps. The method is designed to deal with the 
problem of giving finite and renormalized field 
theoretical definitions of the topological charge 
density, q(x) and of the topological susceptibil- 
ity, x, in lattice gauge theories fl3lfl . The method 
is not new, but is has been recently tested with 
success in certain exactly solvable a-models j}2) 
and, more interesting, it has revealed an unex- 
pected vanishing of \ across the YM deconfining 
temperature |33|] . 

The topological susceptibility has been the ob- 
ject of many theoretical and numerical investi- 
gations not only because it is directly related to 
the mass of the flavor singlet pseudoscalar meson 
(m 2 s — ^prx) H^J], but also because it is a natu- 
ral question to ask whether and how it is possible 
to give on a lattice a notion which would reduce 
to the standard notion of topology in the contin- 
uum limit. For lack of space (see however |35| and 
references therein), I will not mention definitions 
of the topological charge, based on purely geo- 
metrical considerations (the first of which dates 
back to [p6[), nor the construction of Xj based 
on the use of the flavor singlet axial WT1 pTj. 



These two topics are, anyway, outside the aim of 
this talk, as in both cases, by construction, no 
renormalization is expected to be necessary. 

The formula for the topological charge in a con- 
tinuum gauge theory (Pontrjagin number) is Q = 

2 

J d 4 xq c (x), where q c (x) = ■^2^ i i V pJ^{F ilu F pa ) 
is the topological charge density. On a lattice 
one can use the definition 



q(x) 



1 



±4 

2 4 32tt 2 ^ 

(lupa— ±1 



£/i!vpo-Tr(i :, |tll/ fp (T ) 



where -P M „ is a discretization of F^ v p8[ . It is, in 
fact, immediately seen that, in the naive a — > 
limit, q(x) — > a 4 q c (x) + 0(a 6 ). This means that 
in the continuum (field-theoretical) limit, g% — > 
with <7g ~ 1/loga, we will have 

Z Q (gl)a\gl) qc {x) + 0{a & ) 



q(x) 



(17) 



From the previous equations it naturally follows 
that the lattice version of the continuum topo- 
logical susceptibility, Xc = J d 4 x(T(q c (x)q c (0))), 
can be taken to be 



X 



(T(q(x)q(0))) 



According to the general rules of field theory, x 
and x wu l be related in the continuum limit by 

X = Z 2 Q {gl)a\gl)xc + C x (g 2 ) + 0(a 6 ) 
2 



(18) 



The factor Z^(g 2 ) comes from (U7j), while the 
subtraction term, C x (g$) : is a consequence of the 
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mixing of x with the operator Tr (i^i 7 ^) and the 
identity. A NP estimate of Zq^q) and C x (g1) 
relies on the simple observation that short range 
fluctuations (at the scale of the cut-off) are re- 
sponsible for renormalization effects, while phys- 
ical effects, like confinement, come from much 
larger distances (of the order of the correlation 
length, £). When approaching the continuum 
limit, the two scales are widely separated. Using 
a local Monte Carlo updating algorithm, fluctua- 
tions at distances ~ a will be soon thermalized, 
whereas fluctuations at the scale of £ are critically 
slowed down. For a standard local algorithm, like 
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the initial value Q = 1 is reobtained, thus show- 
ing that the global properties of the successive 
gauge configurations one has gone through were 
left unchanged by the heating process. Once Zq 
is known, C x is determined starting from a (topo- 
logically) trivial gauge configuration and measur- 
ing x for a few Monte Carlo steps, until it reaches 
a plateau (Fig. 0). Since in ( |l8| ) the first term is 
zero (Q = 0), the height of the plateau is pre- 
cisely C x . Cooling back the gauge configuration, 
one checks that no unwanted non-trivial contribu- 
tions to x have been introduced during the heat- 
ing process. 

3 



10 20 30 40 

n 




Figure 6. Q as a function of the heating step, n. 



Figure 7. x as a function of the heating step, n. 



Metropolis, the time (# of sweeps) necessary to 
thcrmalize fluctuations at distances d grow pro- 
portionally to d z with z ~ 2. Changes in the 
global (topological) properties of the gauge con- 
figuration are expected to require a much longer 
(exponential) time. To measure Zq a gauge 
configuration with topological charge 1 is placed 
on the lattice and, at any given temperature, 
(3 — 2N c /gg, it is progressively heated (only link- 
variable changes that make the action increase 
are accepted) for a few Monte Carlo sweeps. Zq 
is the value at which Q, measured as a function 
of the Monte Carlo step, n, reaches a plateau 
(Fig. H). Data have very small errors and Zq 
can be extracted with remarkable accuracy. The 
last points of Fig. || (shaded area) are obtained 
by cooling back the configuration. It is seen that 



6. Short distance expansions, higher twists 
and condensates 

Short Distance Expansions (SDE's) are among 
the very few approaches, besides instanton calcu- 
lus and SPT's, that may provide pathways to NP 
analytical calculations in (non-supersymmetric) 
QCD. It is then of the utmost importance to un- 
derstand to what extent corrections to the leading 
term of the expansion can be rigorously defined 
and reliably computed. The heart of the prob- 
lem is that, while the calculation of Wilson co- 
efficients is necessarily truncated to some finite 
order in a s , exponentially small (~ exp(— c/a s )) 
non-perturbative corrections (as those given by 
next order terms) are retained. 
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The whole issue is complicated by the diffi- 
culties arising from the existence of renormalon 
ambiguities in the Borel resummation of PT. As 
is well known [ ^9[ , perturbative series in field 
theories are only asymptotic and in the most 
interesting cases (gauge theories) are not Borel 
summable, as a consequence of the presence of 
singularities located on the positive real axis of 
the Borel w-plane. We recall that a pole at the 
position u = uq > leads to an ambiguity of the 
kind exp(— -%u ) in the Borel resummed series. 

In all instances where a large scale expansion 
is performed (e + e _ -annihilation, DIS scattering, 
Heavy Quark Effective Theory (HQET),...), fur- 
ther renormalon singularities appear. As in the 
full theory, they come from the factorial growth of 
the weight of certain classes of diagrams when the 
running coupling g 2 (k) (k is the loop momentum) 
is expanded as a power series in g 2 (fJ.) (/i is some 
fixed renormalization scale) ^J. The situation can 
be summarized as follows. 

• Borel poles coming from the loop integration 
region where — 6olog(fc//i) > 0, i.e. where k < ji 
(bo > in our conventions), are called Infra- Red 
(IR) renormalons and affect Wilson coefficients. 

• Poles coming from the loop integration re- 
gion where — b \og(k/fi) < 0, i.e. where k > fi, 
are called Ultra- Violet (UV) renormalons. They 
show up in the hadronic matrix elements of the 
Wilson operators. 

A lot of interesting work has been recently 
done (for a nice review see Q) on the ques- 
tion of the appearance and cancellation of renor- 
malon ambiguities. By matching the Wilson ex- 
pansion against full theory calculations, it was 
shown that SDE-renormalon ambiguities must 
cancel between Wilson coefficients and hadronic 
matrix elements. The way it happens is rather 
involved and the details depend on whether one 
uses a soft (e.g. dimensional) or a hard (e.g. lat- 
tice) regularization. 

In spite of this nice result, it has been argued 
in pit] that a further strong assumption on the 



6 There exist also Borel singularities coming from the fac- 
torial growth of the number of diagrams with the order of 
PT. They lie on the negative side of the real axis (think of 
the typical case of \<f> 4 ) and do not affect Borel summabil- 
ity properties. They are not of interest in this discussion. 



relative magnitude of perturbative corrections vs 
NP terms is necessary and unavoidable, if one 
wants to give a precise meaning to large scale ex- 
pansions beyond the leading term. To illustrate 
the point I will closely follow ref. 41 1. 

Let T fl (Q 2 ) = (f\P(Q 2 )\i) be the Fourier 
transform of some bilocal T-product of operators. 
The expansion of Vfi(Q 2 ) for large Q 2 has the 
form 



VftiQ 2 ) = Ci(%){f\Oi(p)\i)- 



C a (Sfr)</|O a (M)|t) 



(19) 



where 0\(ll) and 02(a*) are local operators nor- 
malized at scale with dim02 = dimOi + In. 
For simplicity only two terms have been included 
in ([H]), but the argument can be readily extended 
to any number of terms. The operator P is renor- 
malized (in the full theory) at a scale M and the 
dependence of C\ and Ci on M is understood. 
The usefulness of the expansion ([l9|) lies in the 
fact that short distance NP effects are contained 
in the matrix elements of the operators Oi , while 
the Ci can be reliably computed in PT at large 
\x. It is assumed that the matrix element of 0\ is 
exactly known, either because it is the identity or 
because it is a conserved operator. The question 
is whether one can unambiguously define O2 by 
matching the large Q 2 (perturbative) computa- 
tion of the rhs of (|l9|) in the full theory with the 
form of the expansion in the lhs. Barring excep- 
tional cases, the interesting situation is when O2 
mixes with 0\. 

If dimensional regularization is used, IR renor- 
malon singularities will appear in the Borel trans- 
form of C\, making its perturbative series non- 
Borcl summable. Matching implies that a com- 
pensating UV renormalon ambiguity will affect 
the matrix elements of 02- This simple argu- 
ment already shows where the heart of the prob- 
lem with the definition of higher twist operators 
(such as O2) lies: the perturbative series for the 
Wilson coefficients have renormalon ambiguities 
that are of the same (exponentially small) mag- 
nitude of the NP effects described by the higher 
dimensional terms one is including in the Wilson 
expansion. 

One might imagine to overcome this difficulty 
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by alternatively trying to define higher order 
terms in SDE's by either comparing different ex- 
pansions involving the same higher dimensional 
operators (at the expenses of a certain number of 
predictions) or by resorting to a NP, say lattice, 
determination of the matrix elements of the rele- 
vant higher twist operators. It turns out Q that 
the two methods are equivalent and offer only a 
partial solution of the problem, in the sense that 
1) - after renormalizing O2, C\ is free of the 
leading renormalon ambiguity of order j^yr ~ 



4tt 



r2n) 



2) - Extra renormalon poles, located further 
away along the positive real Borel axis are as- 
sumed to be related to the presence of even higher 
dimensional operators in (^) and could be in 
turn eliminated by sacrificing extra physical pre- 
dictions to fix their ambiguous matrix elements. 

3) - However, since the Wilson coefficients can 
only be known up to a certain order in PT (say, 
up to order (a s ) k ) and the cancellation of the 
leading renormalon starts to be numerically ef- 
fective at some large order in PT (the higher the 
dimension of the involved operators, the larger 
the order the cancellations start to be operative) , 
we must require higher order perturbative contri- 
butions to be negligible with respect to the ex- 
ponentially small terms we are retaining in the 
expansion. 

In formulae this means that for the expan- 
sion (|l9|) to be useful the inequality 



exp[ 



47r I dimo 2 -dimoih 

&oa 3 (Q 2 ) J 



(20) 



must hold [J Unfortunately, whether or not in- 
equalities like (20) are numerically satisfied in real 
life cannot be decided a priori. 

The situation is even more troublesome in the 
case of "condensates", i.e. when one is dealing 
with the vacuum expectation value (v.e.v.) of 
the T-product of, say, two currents, like in e + e~- 
annihilation. The v.e.v's of the local operators 
appearing in the SDE are the NP quantities one 
would like to define properly and utilize to make 



7 Technically the cancellation of the leading renormalon 
ambiguity in Ci ensures that the factorially growing per- 
turbative tail one is neglecting does not sum up to an 
exponential factor larger than the one appearing in (EOh. 



predictions in other physical processes. Naively, 
i.e. according to what we were taught in our grad- 
uate courses in field theory, when a local opera- 
tor, O, can mix with the identity, its v.e.v. will 
be power divergent and should be subtracted out. 
Furthermore, whatever the subtraction prescrip- 
tion is, physical results should not depend on it. 
Thus without loss of_generality, one can always 
use the prescription = — (0)1. The conse- 
quence of this argument is that the v.e.v. of a 
local operator cannot have a physical meaning. 
The only exception is when the v.e.v. in question 
plays the role of order parameter of some symme- 
try. The most obvious example of this situation 
is (tpip), which is the order parameter of chiral 
symmetry. In this case one can use the WTI of 
the chiral symmetry to define what is to be meant 
by (ipip) and, if required, to relate its value to the 
mass of the Goldstone boson, when m q ^ 0. 

This situation should be contrasted with the 
case of the condensate F 2 = (Tr(F M „F M „)). F 2 
is not the order parameter of any symmetry and, 
in fact, it does not appear in any useful WTI. 
Besides, if we recall that the v.e.v. of the trace 
of the energy-momentum density tensor, (9^), is 
proportional to F 2 , we would be very strongly 
tempted to conclude that F 2 = 0, to prevent 
the vacuum to have infinite total energy. The 
last statement follows from the observation that 
Lorentz invariance implies (0^) — kg^, so that 
vacuum energy finiteness requires (#00) = 0, lead- 
ing to the conclusion k = and hence F 2 = 0. 
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